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Abstract In this paper, limitations of the common method measuring intrinsic spatial resolution of the GEM 
imaging detector are presented. Through theoretical analysis and experimental verification, we have improved 
the common method to avoid these limitations. Using these improved methods, more precise measurement of 
intrinsic spatial resolution are obtained. 
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1 Introduction 

The GEM(Gas Electron Multiplier) is a typical 
Micro-Pattern gaseous detector, first invented in high 
energy physics, then applied in many other fields. We 
have constructed a 2-D imaging detector using triple- 
GEM for BRSF(Beijing Synchrotron Radiation Facil- 
ity) with active area of 200inmx200mm, which is able 
to detect X-ray with high spatial resolution(Fig. [T]). 
The GEM consists of a thin, metal-coated polymer 
foil, etched by high density of holes. When potential 
is applied between electrodes up and down, radiation 
electrons drift into the holes, multiply and transfer to 
the other side. Each hole acts as an individual pro- 
portional amplifier. In our detector, three GEM foils 
are used to get much higher gain. 

One important specification of an imaging detec- 
tor is the intrinsic spatial resolution(resolution will be 
used for short later in this paper). There are several 
definitions of resolution according to different crite- 
ria, however, all of them are equivalent in fact. Take 
FWHM and the standard deviation (Tq for instance: 
when the consideration function is the normal distri- 
bution, the relationship between standard deviation 



(To and FWHM is FWHM = 2y/2ln2ao. In this pa- 
per, the standard deviation is chosen to be used as 
the resolution [l[ Q. Overall, resolution is to describe 
the resolving power of imaging detectors, so how to 
measure the resolution accurately is important. 
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Fig. 1. The schematic view of our GEM detec- 
tor structure. 
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2 Analysis of the measurement of spa- 
tial resolution of imaging detectors 

2.1 The common method to measure spatial 
resolution 

The common method to measure the resolution is 
usuahy to let beams collimated by a slit (sometimes 
by a hole or a blade) (Fig. [2]). The measurement result 
is always a distribution, of which standard deviation 
is a. Then the resolution of the detector, of which 
standard deviation is Uo, is obtained by Eq. ([1]). 



(7o = V^^^ (1) 

where h is the width of the slit. 




Fig. 2. The schematic view of beams collimated 
by a slit. The X-ray beams go through the 
collimator, which is horizontal. The detector 
surface is vertical. 



We have done the measurement using this method 
in the condition: d — 20mm, h — 0.2mm, and 
I = 40mm. The a of experimental data distribution 
is 1.33mm. The resolution is calculated by Eq. (HI). 

do = Vl-33' - 0.22 = i.3i49TOm (2) 

Through a careful study, we have found limita- 
tions of this method and the error introduced in this 
methodfSec. 123)) . 

2.2 Resolution of the detector, beam inten- 
sity distribution and experimental data 
distribution 

When we use the GEM detector for X-ray imaging 
of an object, the beam(after going through the colli- 
mator) intensity on surface of the detector is always a 
random variable, of which distribution is called beam 
intensity distribution. The distribution of data 
measured by the detector is called experimental 



data distribution. Due to the electrons diffusion 
in the detector and error existing in the measure- 
ment, the experimental data distribution can not re- 
flect the resolution function accurately. Experimental 
data distributions are often the result of the resolu- 
tion function that is modified by the beam intensity 
distribution. In fact, the experimental data random 
variable Z is a superposition of the beam intensity 
random variable X and detector resolution random 
variable Y as shown in Eq. ([3]). 

Z = X + Y (3) 

According to the central limit theorem, the res- 
olution function usually follows the normal distribu- 
tion. Hence, in the most general case, experimental 
data distributions are described by a convolution of 
a beam intensity distribution and the detector reso- 
lution function, where the p. d.f. (probability density 
function) of experimental data distributions is g{x'), 
that of beam intensity distribution f{x) and that of 
resolution function r{x' — x) (Eq. (|4])). 

g{x')^ I r{x' -x)f{x)dx (4) 

fix is the domain of the beam intensity random vari- 
able X. Its specific form is shown in Eq. ^ and 
Eq. ©. 

Hence, the resolution can be calculated by fit- 
ting the experimental data distribution with g{x') 
(Sec. 13. 1[) or calculating r{x' — x) by deconvolu- 
tion(Sec.[3l2l). 

2.3 Limitations of the common method 

With different beam intensity distributions, differ- 
ent experimental data distributions are obtained ac- 
cordingly, as shown in Table [TJ From the table, only 
when the beam intensity distribution and the detec- 
tor resolution function both follow the normal distri- 
bution, can the resolution be calculated by Eq. ([1]). 
In many situations the beam intensity distribution 
does not follow normal distribution. In these condi- 
tions, if the width of the slit is smaller than tenth 
of the resolution(FWHM), the common method is 
still apphcable according to GB/T 18989-20cl3. This 
is the scope of application of the common method. 
To avoid this limitation, methods are improved to 
get more precise resolution. Experiments with GEM 
detectors 0, 0| have been done to confirm these im- 
proved methods(Sec. [Sj. 
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Table 1. Experimental data distributions as convolution of resolution function and different beam intensity 
distributions. 



p.d.f. of resolution function 



p.d.f. of beam intensity distribution 



p.d.f. of experimental data distribution 



/27r(7o 



(x'-i)2 



f{x) = 5{x-xo) 
(extremely narrow slit) 

(uniform distribution; a, b is the lower and upper limit) 



9(x-') = 



■\/27r(To 



e 2<7o2 



-e 



(normal distrubution) 



gix') -. 



V2ir(h2+CTo2) 



e 2(fi2+„^2) 



* -F(a;) is cumulative distribution function of f{x). 



3 Improved methods to measure the 
spatial resolution of imaging detec- 
tors 

3.1 Improved method 1: using a slit as the 
collimator with convolution fit 

As shown in Fig. [21 the surface of detector, which 
is d + / away from the beam source, is perpendicular 
to the beam line. So Fig. [2] actually shows the one 
dimension projection of the measurement frame. The 
beam source is a narrow strip source limited by the 
slit with a definite emission angle associated with the 
total domain il^ in Eq. (|4]). The total effect is that 
the intensity along the X-axis on the surface of the de- 
tector is proportional to the range of the source from 



where the X-ray beam can reach the surface. Due to 
the effect of the collimator, the beam intensity distri- 
bution is divided into three parts. P.d.f. of the beam 
intensity distribution is as shown in Eq. ([5]). 

h+fix-fio + 'l) xe[-{^ + h^)+fio, 

(5) 

where /io is the coordinate of the middle of the slit. 
f{x) should be normalized. 



r{t+h^)+Mo 

-(^+hi) + MQ 



fix)dx-- 



dx+ 



r+Mo 



h dxH 



r+Mo 



(| + hi) + MO 



d h 
h+-[x-Ho+^) 



dxj =1 
(6) 



From Eq. ([6]), we can get the normalization coef- 
ficient. 

a^l/h\-^ + l) (7) 

It is worth to remind that if without the slit, emission 
angle of the line source is 2ti. 

Curve of the beam intensity distribution, solid 
black line with shape of trapezoid (like a dam), is 
shown in Fig. |3l When CTq is very small, the experi- 
mental data distribution is close to the beam intensity 
distribution(dashed line in Fig. [3|). Otherwise, when 
(To is very large, experimental data distribution do 
not reflect the beam intensity distribution but just re- 
flect detector resolution function r{x' —x) itself (dash- 
dotted line in Fig. [3]). 
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Fig. 3. Simulation study: the experimental 
data distributions as a function of the same 
collimator but different resolution of the de- 
tector. 
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In other words, if the detector resolution is much 
larger than the slit width h, g{x') is effectively ap- 
proximated by r(x' — x) and the collimator effect 
can be ignored, so an extremely narrow slit can be 
used as a collimator to verify these improved meth- 
ods(Sec.[4l). 



X-ray Source 
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Readout 
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Fig. 4. The schematic view of the incident an- 



gle affection. 

Attention should be paid to one more affecting 
factor of oblique incidence shown in Fig. 01 Because 
of the incident angle Z6', the quantity of the electrons 
on the top of GEMl follows the uniformity distribu- 
tion over the domain of a tan (a is the space of the 
drift area of the GEM detector) instead of a point. 
So, all of the measurement results of resolution in 
this paper should remove the incident angle affection 
by subtracting atan^. 

In this method, the beam intensity distribution 
is known, the resolution function follows the normal 
distribution, of which the standard deviation CTq can 
be treated as an undetermined coefficient. When we 
measure the resolution in experiment, the convolution 
function(Eq. ([8])) is used to fit the experimental data 
distribution [il to obtain the resolution of the detec- 
tor. The measurement result on GEM detector using 
this method is determined to be ctq = 65.0/xm. 



+;4)+Mo 1 



e 2,^0^ 



r+Mo 



f 2+ ^ 



(x'-x)2 

2-o2 -h Ax 



| + /io,| + Mo) 



f+MO 



-(f + h^) + MO V27r<To 



[/i + f(x-Mo + t)]cla: xe [-(f + /ii) + Aio,-HMo] 



The normalization coefficient q; = 1//i^(^ + 1) 



(8) 



3.2 Improved method 2: using a blade as the 
coUimator with deconvolution and convo- 
lution fit 

In a 2-D image obtained from an imaging detec- 
tor, it is reasonable that the sharper the image edge 
is, the more precipitous its projection histogram is. 
As shown in Fig. [SJ the beam intensity distribution 
of the edge is a step distribution if a blade is used to 
cut the beam. The beam intensity above the blade 
nearly follows the uniform distribution. The beam 
intensity below the blade is zero. As the resolution 
function is a Gaussian distribution, the experimental 
distribution is superposition of the step distribution 
and Gaussian distribution, which is in fact the cumu- 
lative Gaussian distribution!] 



g(x')= / r{t)dt (9) 
The resolution function can be obtained by deriva- 



tion of X on Eq. (j9]). do of the cumulative Gaussian 
distribution is the resolution. The process of solv- 
ing the derivation(gradient) is that of deconvolution 
too. The measurement frame on GEM detector us- 
ing this method is shown in Fig. [5] and the result is 
tJo = 71.3^TO(Fig. [6|). It is a little larger than the re- 
sult of the first improved method, because the beam 
lines are not absolutely parallel. 




detector 



Fig. 5. The schematic view of using blade as 
the collimator. d = 2cm, / = 4cm. 
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Experimental data distribution 



Derivation of experimental distribution tiistogi 
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Fig. 6. Measurement result, (a) is the experi- 
mental data distribution; (b) is the derivation 



of the experimental data distribution. 

It is worthy to note that this method requires large 
statistics, otherwise, if there is no sufficient statistics 
per bin, taking the Gaussian distribution for an exam- 
ple, the derivation distribution of the same Gaussian 
p.d.f for different bin width will be different, as shown 
in Fig. [71 which well yield different (Tq for Gaussian 
p.d.f. In our case, the total statistic is 2 x 10^, while 
the bin width is lOfim. 



Histogram of Gauss distribution 




Derivation of Gauss liislogra 




with bin width 0.12 Derivation of Gauss histogram with bin width 0.03 




Derivation of Gauss histogram with bin width 0.006 




(a) X (b) X (c) X (d) X 

Fig. 7. Limited statistics in a bin of a histogram for a Gaussian distribution will yield different derivation 
distributions, (a) is a Gaussian distribution; (b), (c), (d) are derivation of the Gaussian distribution with 
different bin-width. 



To avoid this kind of drawback, we improved this 
method, measurement architecture of which is the 
same to the deconvolution method but with a dif- 
ferent data processing. As described above, the edge 
of an image follows cumulative Gaussian distribution. 
Indeed, the standard deviation of this Gaussian distri- 
bution is the resolution ctq. So the cumulative Gaus- 
sian distribution with CTq as the undetermined coeffi- 
cient can be used to fit the experimental data distri- 
bution. In this way, the resolution is determined to 
be (To = (Fig. |8]). 
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Fig. 8. The experimental distribution using a 
blade as the collimator. 



Verification the improved methods 
by using a extremely narrow slit as 
the collimator 



As mentioned in Sec. 12. 3[ if the width of the 
collimator is narrow enough, less than tenth of the 
resolution(FWHM), the common method is applica- 
ble. In this method, an extremely narrow slit is used 
as the collimator. The measurement architecture is 
shown in Fig. [2l where h = 0.01mm, I = 40mm and 
d = 20mm. Because the slit is so narrow that can it 
allow a very little part of beams to go through. High- 
intensity beams are required to get sufficient statis- 
tics. It is almost impossible to use the X-ray tube as 
the beam source in this situation. The measurement 
has been done at BSRF which can provide enough 
high intensity beams. 
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Fig. 9. The measurement result using ex- 
tremely narrow slit as the collimator. 
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The result is shown in Fig. |9l In this way, the res- 
olution is (To — 59.9/im. The result strongly validates 
the consistence with those of the improved methods. 

5 Measurement using the Rayleigh 
criterion 

The Rayleigh criterion is the generally accepted 
criterion for the minimum resolvable detail when the 
first diffraction minimum of the image of one source 
point coincides with the maximum of another. That 
is to say, the intensity of the saddle between two peaks 
of the points is 81% the intensity of each peak(3]. 



Experimental data distribution Resolution and beam intensity distribution 




Fig. 10. The measurement result using the col- 
limator with 3 slits. Left: the measurement 
data distribution; Right: the beam intensity 
distribution and the resolution function. 



To use the Rayleigh criterion, we have measured 
the distribution using a collimator, which has 3 slits 
with 0.32mm in pitch. Each slit is 0.20m7Ti wide. The 
measurement result is shown in Fig. 1101 As the width 
of each slit can not be ignored, the minimum resolv- 
able pitch is 0.32mm(Fig. [TUl left). After removing 
the effects of width of each slit by deconvolution fit, 
the resolution is ao — 76.2/im. It is a little larger than 
the results of the first and second improved methods 
mentioned above. That is because of the definition 



of the Rayleigh criterion. If the considered function 
is the normal distribution, according to the Rayleigh 
criterion, when the intensity of the saddle between 
two peaks of the points is 81% the intensity of each 
peak, the distance center to center of the two peaks 
is 2.69(To. Yet, the distance of the first and second 
improved methods is the FWHM, which is 2.355cro. 
So the resolution using the Rayleigh criterion is 1.142 
(1.142 = HI?) times larger than those of the first and 

^ 2.355 ^ ^ 

second improved methods. 

6 Summary and conclusion 

The imaging detector's resolution can be limited 
by diffraction causing blurring of the image. By care- 
ful study of the relationship between the spatial reso- 
lution and the collimator, we give the scope of ap- 
plication of the common method and improve the 
method. The measurement using the improved meth- 
ods gives more precise spatial resolution. Experimen- 
tal validation has been done on the GEM detector as 
shown in Tab[2l 

Table 2. Summary of results of the measuring 
methods. 



Method# (TO with incident angle correction 



Common method 


1314.9 nm 


Method 1 


65.0 fim 


Method 2-1 


71.3 fim 


Method 2-2 


63.3 fj,m 


Extremely narrow slit 


59.9 fim 


Using the Rayleigh Criterion 


66.7 fim* 



* 66.7 = 76.2/1.142. 



The measurement is an exploration of the mea- 
surement of spatial resolution and a general reference 
for that of other imaging detectors. Further study 
taking more factors in to account is in progress. 
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